A square root of the harmonic oscillator 



R Delbourgo* 

Physics Department, University of Tasmania, 
GPO Box 252C Hobart, Australia 7001 
(February 1, 2008) 

Abstract 

Allowing for the inclusion of the parity operator, it is possible to construct a 
model of an oscillator whose Hamiltonian admits an exact square root which 
is rather different from the conventional approach based on creation and an- 
nihilation operators. We outline such a model, the method of solution and 

some generalisations. 
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*In memory of Lorella M Jones who sadly passed away in February 1995 and to whom the original 
idea in this paper is due. 

1 



I. INTRODUCTION 



The simplest way of solving the harmonic oscillator uses the fruitful concept of annihila- 
tion A and creation operators. A and A* are of course non-hermitian and therefore differ 
from one another just like a complex number and its conjugate; nonetheless, in a loose way, 
they can be construed as providing a square root of the harmonic Hamiltonian, since the 
number operator is given by N = A^A. The same idea applies to supersymmetric models 
where the Hamiltonian is expressed as a quadratic function of hermitian fermionic operators 
tj, viz. V, ; r'' /•;/•• . In this paper we wish to exhibit a model for which there exists an exact 
square root of the Hamiltonian that is not far removed from the one for an oscillator. In that 
sense it represents a much closer analogy with the Dirac equation, which is a square root 
of the Klein-Gordon equation at the price of doubling the number of components (namely 
introducing spin). Indeed, we shall find something similar in the system described below. 
To dispose of unwanted constants and thereby simplify the algebra, let us set the mass and 
circular frequency u of the oscillator equal to 1, as well as h. 

It is worth mentioning that a number of papers have appeared concerning the Dirac 
oscillator. These correspond to making a minimal substitution p — > p — if3r'm the free Dirac 
equation and they can be regarded as placing the spinor particle in a linearly increasing 
potential. We should stress that our approach is quite different from those papers since it 
even applies to a bosonic system. The crucial respect where we differ is our reliance upon 
the parity operator Q. 

The basic point is to notice that whereas A = P + iX is non-hermitian, the combination 
D = P + iQX is real ]|]. Here P stands for the momentum operator, X is the position op- 
erator and we are working in one spatial dimension for the moment. Using the fundamental 
commutators, 

[X,P]=i, {X,Q} = {P,Q} = 0, Q 2 = l, 
and squaring D, one gets 

H = D 2 = P 2 + X 2 - Q = 2H osc - Q. (1) 
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Thus apart from a parity term Q ||, which commutes with H osc , we get an exact square root 
of the oscillator. It is an unexpected and striking result that we do not require non-hermitian 
operators for 'square-rooting' model (1). 

In the next three sections, we shall study the eigenvalue equation for D, extend the 
idea to higher dimensions and generalise the analysis to the Dirac equation. A full set of 
physical solutions will be determined which do not depend on an external vector; we shall 
also investigate zero eigenvalues that are characterised by some external (null) wavenumber. 

II. EIGENFUNCTIONS OF D 

Because the two parts of H behave differently under parity, we may anticipate that a 
doubling of components is needed and that is precisely what the following representation of 
the various operators in H gives, 

P — > — io\ d/dx, X — > aix, Q — > cr 2 - (2) 

Thus D is represented by the matrix, 

D —— > — io\ d/dx + (J3X 



( x —id/dx 



—id/dx —x j 

Letting (f(x),g(x)) denote the eigenfunctions of D, it pays to form the combinations f± 
J ± ig, corresponding to projections (1 =F Q), so as to obtain the two coupled first order 
equations, 

(±x + d/dx)f T = ±A/±, (3) 

where A connotes an eigenvalue of D. In turn these can be converted into uncoupled second 
order equations, 

f±"-(x 2 -\ 2 ±l)f± = 0. (4) 
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The standard method of solution, which is to put f± = F± exp(— x 2 /2) and work out the 
(terminating) series expansion of F±, shows that the excited states are 

/+ oc H N exp(-x 2 /2), /_ oc ±H N+1 exp(-x 2 /2); 



with A = ±y 2(iV + 1), N = 0, 1, intrinsic parity Q = (— 1) N+1 ; in contrast, the ground 

state is 

/ = -i g oc exp(-x 2 /2); A = 0, 

with intrinsic parity Q — 1. Here H N {x) signify the usual Hermite functions. Altogether 
then, with proper normalization, we write the (doubled) excited states in the form 



ip N (x) 



-texp(-x 2 /2)_ ( i(H N -i ± H N /V2N) ^ 
K (H N ^ t H N /V2N) } 



(5) 



with A = ±\ / 2N; N = 1, 2, . . .. The ground state (A = 0) is 

(A 



exp(— x 2 /2) 
^27? 



(6) 



As a matter of fact one may regard the ground state as the continuation of (5) to N = 0, 
remembering that (—1)! = oo. 

One can then go on to determine expectation values, uncertainties and such like, using 
well-known properties of Hermite functions. For instance, the ground state quite obviously 
yields, 

(X) = (P) = 0, AX.AP = 1/2, 
whereas the excited states (N — 1, 2, . . .) have 

(X) = 0, (P) = ±^/N/2, AX.AP = N/V2, etc. 

However instead of elaborating upon this 1-dimensional model, it is much more interesting 
to examine its higher-dimensional generalizations. 
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III. HIGHER DIMENSIONS 



With several independent coordinates (j = 1, . . .,D) we can contemplate at least two 
distinct choices for the square roots Dj, viz. Pj + iQjXj or Pj + iQXj, where Qj reflects Xj 
alone, whereas Q = T7j Qj is the full inversion operator. The problem with the first choice 
is quite serious: the components Dj do not transform as a true vector under rotation, in 
contrast to the second choice. However in the latter case one should observe that 

[Dj, D k ] = [Pj + iQXj, P k + iQX k ] 

= 2iQ(X j P k - X k Pj) = 2iQL jk . (7) 

Nevertheless it remains true that J2 D 2 is effectively the D-dimensional oscillator, barring a 
parity term: 

DjDj = PjPj + XjXj - DQ = 2H osc - DQ. 

We will therefore stick with the second choice. 

In order to determine the eigenfunctions of this Hamiltonian, we cannot take simultane- 
ous eigenfunctions of the Dj separately since the components will not commute. Still, we 
may construct joint eigenfunctions of DjDj, Q and L 2 , where L jk are the D(D — l)/2 com- 
ponents of the angular momentum tensor. Now in D-dimensions, proceeding to appropriate 
spherical coordinates, 

^{x) = R(r)Y(x), Lj k L jk .Y{x) = 21(1 + D- 2).Y(x), 

the radial equation for D 2 i\) = \ 2 ^ can be reduced to O r R(r) = 0, where 

d 2 D-ld e(£ + d-2) 2 x2 ^ 
O r = — + ^ '- —r+X+ DQ. 8 

dr z r dr r z 

As ever we look for a solution, R(r) = H(r) exp(— r 2 /2) where 7i is a terminating power 
series in r and obeys V r H(r) = 0, with 

^ d 2 , D ~ l , 2 „ M + D-2) 
D r = -— + 2r — + A + DQ — D — -± 

dr A r dr r A 
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It is an elementary exercise in differential equations to prove that 1~i{r) = Y^j c j r with 
c 7^ where the coefficients satisfy the recurrence relation, 

(j + 2)(j + 2£ + D)c j+2 = (2j + 21 + D - DQ - X 2 )c r 

Clearly we must set the odd terms, starting with ci, equal to and require that 

A 2 = 4n r + 2£ + D - DQ; (Q = ±1), 

where the radial quantum number n r = 0, 1, 2, . . . has the same range of integer values as I. 
In short we have found the complete solution to this problem as a 'radial Hermite' function 
multiplying a spherical harmonic and radial gaussian, whose coefficients are fixed by 

c j+ 2 _ 2(j - 2n r ) 
c 3 ~ (j + 2)(j + 2£ + dy 

Finally, the eigenvalues A 2 of D 2 are given by 

X 2 = 2N = A(2n r + £) + (D t D); Q = ±1. (9) 

It is tempting to suppose that the above discussion applies to the relativistic Klein- 
Gordon operator in the particular case where D — > 4 and we continue in the time coordinate, 
x 4 — > it, Pi — > iE. Although it is true that the solutions 

<Mx) = H(p).Y e (p) exp( V/2); P 2 = -x 2 , (10) 

with 

(A 

A 2 =4n r + 2£+ , for Q = ±1, 

are certainly normalizable over 3-space at t — 0, they have unacceptable behaviour in the 
distant past or future and share some of the time-dependence problems of Bethe-Salpeter 
wavefunctions that are not evaluated in the instantaneous approximation. We therefore 
somewhat hesitant to claim || that such a covariant Hamiltonian (= D^D 11 ) with its Gaus- 
sian factor exp(x 2 /2) makes good physical sense. 
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IV. ANOTHER DIRAC OSCILLATOR 



We are now going to extend the discussion covariantly to the Dirac equation, but differ- 
ently from refs |l]-f|]. Continuing to define Q as the total inversion (x — > —x) operator, the 
Dirac oscillator equation 



1 .(P + iQXW = i\4>, 



(11) 



can be represented in the doubled 4-component form 7.(<9<Ti + ixa^jif) = or 



7.x i'j.d 



\fA (A 

= i\ 



f 
\ 9 J 



(12) 



It helps to form the combinations f± = / ± ig as before, so that the equations simplify to 
the coupled pair 



j.(d±x)f T = =FiA/±. 



(13) 



"Squaring" , we arrive at the uncoupled pair of differential equations, 



(d — x ± 4 ± a^ v L^ u ) f T = X f_ 



(14) 



since jk ty- fi/c = 7 / V(x M <9 i , - d^x u ) = -4 - cr^L^ with = i(Xfj,d u - x u d^). 

Now products of spinors and 4-dimensional spherical harmonics decompose relative to 
the Lorentz group as the sum of two parity conserving irreducible representations A and B 

4,o)e(o I i)]®(^) = [A]e[s], 



where 



.4 
B 



+ 1 L 
2 ' 2' 



.£ i + 1 



- 1 £ £ e-i . 



The square of the total 4-dimensional angular momentum J = L + a/2 acts thus: 



J^J IMU .[A] = (2f + 6i + 3).[A], 



(15) 



whereupon it follows that 

a^L^.[A] = 2£.[A]; 

similarly for the wavefunction [B], upon effecting the substitution t — > i — 1. Hence on one 
of these irreducible wavefunctions (labelled by A say), /a+(x) = -^a+(p)-^a(p) our oscillator 
equation simplifies to 

(d 2 -x 2 -4-2£-\ 2 )f A+ (x) = 0. (16) 

As usual, put Fa+(p) = 'Hip) exp(— p 2 /2); p 2 = — x 2 , to get a radial equation O p H(p) = 0, 
having 

The (polynomial in p) eigenfunctions are readily found, 

Up 

^4 + (p) = E c ^ 

j=0 



with A 2 = 4n p ; n p = 0, 1, . . . and 

c 7^ 0; j even. 



Cj (j + 2)(j + 2£ + A) 
One may carry out a similar analysis of the equations; this time one ends up with 
the eigenvalues A 2 = An p + 2i + 8. Of course we get exactly the same answers for the [B] 
eigenfunctions, apart from the change t — > £ — 1. This then completes the analysis when 
there is no dependence of /, g on an external wave vector. 

Let us now have a brief look at the zero eigenvalue problem allowing for some external 
wavenumber — which is therefore analogous to seeking plane wave solutions of the ordinary 
free Dirac equation. Recall from (13) that if A = 0, the first order equations decouple. In 
particular the equation for f-(x) has asymptotic behaviour exp(— x 2 )/2 which is certainly 
unacceptable in the space-like direction; that is we are obliged to set /_ = 0. On the other 
hand, putting f+(x) = F + (x) exp(x 2 /2), one simply remains with r y.dF + = 0, corresponding 
to a free massless equation. Hence we arrive at the solution (A = 0), 

f+(x) = [u+ / kvexp(—ik.x)5(k 2 )]exp(x 2 /2). 



In short, the ground state is fixed by a null wave number and may be written 



exp(x 2 /2), 



u+ Jkv exp(— ik.x)5(k 2 ) 
— i[u+ Jkv exp(— ik.x)5{k 2 )} 

where u and v are constant 4-spinors. We have not succeeded in determining solutions for 
A 7^ that are specified in terms of a non-null k, as the equations are much more complicated 
in that case. 

Finally, observe that what progress one has achieved depends critically on the fact that 
there is a double degeneracy when the eigenvalue A is nonzero, but a single degeneracy for 
the ground state A = 0; this happy accident has arisen because the Hamiltonian is exactly 
that of the oscillator plus an integer times the parity. This smacks strongly of the situation in 
super symmetric Hamiltonian models and it would not be surprising if such an interpretation 
can be found. 
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